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ABSTRACT 

We analyse the stability of a magnetized medium consisting of a neutral fluid and a 
fluid of charged particles, coupled to each other through a drag force and exposed to 
differential body forces (for example, as the result of radiation forces on one phase). 
We consider a uniform equilibrium and simple model input physics, but do not arbi- 
trarily restrict the relative orientations of the magnetic field, slip velocity and wave 
vector of the disturbance. We find several instabilities and classify these in terms of 
wave resonances. We briefly apply our results to the structure of SiO maser regions 
appearing in the winds from late-type stars. 
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1 INTRODUCTION 

Multiphase flows are a widespread and important phenomenon in astrophysics. The difference between heating and cooling 
rates in different components of astrophysical gases often leads to the formation of a multicomponent medium, in which several 
phases with widely separate temperatures coexist near to pressure equilibrium. Effective multiphase behaviour can also result 
from differential coupling of distinct particle species to local magnetic fields or radiation driving forces. To first order, these 
differential forces will lead to drift velocities between the different components of the fluid, limited by the effect of frictional 
terms. However, this means that there is local source of free energy in the flow. 

Radiation pressure on dust, for example, plays a major role in many models of the acceleration of winds from highly 



evolved, low-mass stars (e.g. MacGregor & Stencel 1992). The dust streams through the neutral gas and transmits momentum 
to it through collisions. The streaming of dust through neutral gas has also received attention in many other astrophysical 
contexts, including the evolution of dust bounded Hll regions ( [Cochran fc Ostriker 1977 ) and the radiation-driven implosion 
of dense globules ( Sanford, Whitaker fc Klein 198^ ) . Radiation pressure on dust may levitate interstellar clouds above the disc 
of the Milky Way (Franco et al. 1991): in such clouds dust particles will stream through the neutral gas. No magnetic field was 
included in any of these studies. Hartquist fc Havnes (1994) identified conditions under which dust grains are well-coupled 
to the magnetic field when the grains are driven by radiation pressure. In many circumstances the dust, gas phase ions, and 



electrons may be treated as a single fluid. Except for studies of the Wardle instability of shocks in dusty media (Wardle 199C ; 



Stone 1997 



Mac Low & Smith 1997), investigations of instabilities in weakly ionized astrophysical media on lengthscales short 



There is a considerable literature on large-scale instabilities driven by the self-gravity of multifluid media (e.g. 


Mouschovias 


1976; 


Nakano 1976; Huba 1990; 


Brandenburg & Zweibel 1995; 


Balsara 1996; Zweibel 1998; 


Kamaya & Nishi 2000; 


Mamun & 



Shukla 2001 ) . These papers differ in aspects such as the number of different flow components assumed, the precise nature of 
the inter-species coupling and the inclusion of processes such as the self-gravity of the flow and large -scale gradien ts in fl ow 
properties. Many of these papers recover a large scale instability, first described by Mouschovias ( 1976 ) and Nakano ( 1976 ), in 
which the diffusion of magnetic field out of a self-gravitating clump reduces magnetic support, leading eventually to collapse. 
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In some, rapidly growing, small lengthscale instabilities are found (Huba 199C; Kamaya & Nishi 2000; Mamun & Shukla 2001), 



but, to date, rather restricted classes of relative orientation of magnetic field, mean flow and wave vector have been assumed. 

We also note that there are many non-astronomical examples of interspersed multiphase flows, such as clouds, fluidized 
beds and microbial suspensions, which have been studied extensively. For example, Childress & Spiegel (1975) find buoyant 
instabilities, similar to those we discuss here, in systems of finite extent in both astrophysical and terrestrial contexts. 

The past work on the Wardle instability and molecular cloud support has treated inhomogeneous media, in which 
the streaming is induced, for instance, by impulsive acceleration or by large scale variations of the magnetic field. While 
astrophysical flows are necessarily inhomogeneous in the large, these variations can serve to obscure the mechanisms of small- 
scale instability. Given this wide variety of physical mechanism and equilibrium structure, in this paper we outline a general 
analysis for a simplified physical model, in which a charged magnetized fiuid streams through a neutral fluid as a result of 
differential body forces. This model might most directly be related to flows with differential radiative forces on the fluids, but 
can be applied more widely. By assuming uniform initial conditions and treating the modes which we flnd as distributions, we 
can study the stability of short wavelength modes in general, without needing to treat the specific global features which are 
important for longer wavelengths. Our analysis complements the previous work described above, by giving stability criteria 
for wave-vectors of arbitrary orientation and all initial angles between the body forces and magnetic fields, albeit for rather 
simpler input physics. 

In Section ^, we present the basic two-fiuid equations and derive the dispersion relation for linear waves. In Section ^ 
we present numerical solutions of the dispersion relation. For small wavelengths, we find that 'resonances' (where distinct 
modes have similar phase velocities) are important in understanding the stability properties, and discuss a graphical method 
of locating these resonances in general geometries. We then, in Section ^, analyse the stability of the solutions of the dispersion 
relation, proceeding from general analysis to specific analytic stability criteria for short and long wavelengths. These criteria 
compare well with the numerical results in the previous section, and confirm their generality. In Section ^ we apply the long 
wavelength results to the properties of SiO maser spots in late-type stars. Finally, in Section M, we summarize our results. 



2 BASIC EQUATIONS AND DISPERSION RELATION 



In the present paper, we study the stability of two-fluid flows in which one component is coupled to a magnetic fleld. Differential 
forces on the two fluids lead to inter-phase slip in the equilibrium solution. We attempt to characterise the general properties 
which allow instabilities to feed off the slip energy. While the system we consider is simplified, it allows us to analyse the 
processes from which instability results in some detail. 

The equations we treat are those of continuity and momentum for neutrals 

dpn 



^ + V.(p„u„ 



0, 
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and the ideal MHD equations for ions 
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where equation (jH]) maintains the solenoidal condition V.B = so long as it is true initially, and the drag terms correspond to 
Stokes' law with drag coefficient A. Here pn, Un, Cn and g are, respectively, the mean density, velocity, the effective isothermal 
sound speed and the net acceleration for the neutrals; p, u, Ci and gj are the equivalents for the ions; and A is a frictional 
coupling constant. The value of the vacuum permeability p determines the magnetic field unit: conventional values include 
1, 4-K and 47r x lO"'^. These equations are similar to those used in previous work (e.g. Shu 1983; Kamaya & Nishi 2000), 



except that in the present paper we do not restrict the relative orientation of the various vector fields (while, for the present, 
neglecting some of the physical terms included by these earlier authors). 

Two distinct pressure terms are used for the two distinct phases. This might be taken as an assumption that the scattering 
between gas particles of the same phase is far more rapid than that between particles of differing phases. The limits Ci — > 
and Cn — » are relevant in particular contexts, but have degenerate eigenmodes: by assuming finite values of these parameters, 
the degeneracies are lifted within the present analysis. 

The forces on the ionized compon ent are oft e n dominated by th e effects of the overall curvature of the magnetic field, in 
which case one can assume uq.B ~ (3hul98S; Mouschovias 1987 ). It is consistent to study small-scale instabilities in the 
present of such large scale gradients (e.g. Huba 199C ). However, in this case at least part of the equilibrium force on one phase 
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will not be proportional to mass, as we have assumed above: the detailed instability criteria will be somewhat different from 
those which we derive below, but the general behaviour should be similar. In what follows, we study the equations for general 
orientations of the various vector parameters, while bearing in mind the practical importance of cases in which uq.B is small. 

As in most previous papers, we have not included internal viscous terms for the individual phases, although coUisional 
processes will in general lead to viscosity within the phases as well as inter-phase drag. This viscosity will lead to the 
stabilization of unstable wave modes at short wavelengths. We consider the effects of viscous terms in the context of an 
astrophysical example in Section ^, and verify that, in that case, viscosity can be neglected at the wavelengths which interest 
us. We will present a detailed treatment of the stability of flows including both inter-phase drag and internal viscosity in a 
future paper. 

We perturb about an equilibrium with a constant slip velocity uq = u — Un between the phases. The slip velocity satisfies 
the equilibrium condition 

- g + Apuo = -g; - ApnUo = 0. (6) 
Linearizing about this equilibrium with Un = Vn, u = uq + v, p = po + ^, pn = pn,o + 6n and B = Bo + /3, we find 

^n + Pn,0V.Vn = 0, (7) 

Vn = 2_V6ln + Apo(v-Vn)+A6»Uo, (8) 

Pn,0 

e uo.Ve + poV.v = 0, (9) 

1 

v-mo.Vv = -— Ve- A6l„uo + Ap„,o(vn-v) + (V x /?) x Bq, (fO) 

po PPo 

/3 = -uo.V/3-f- (Bo.V)v-Bo(V.v). (11) 
In what follows we will suppress indices on po and Pn,o- Looking for solutions of form expi(k.x — uit), we find 

- Uj9n = -pnk.Vn, (12) 

{-Lu — iXp)vn = ^k6ln — jA(pv -I- U06'), (13) 

Pn 

{uo.k -Lu)e = -pk.v, (14) 

1 

(uo.k - - iAp„)v = — !-ke-iA(pnVn-uoen) + — [(Bo.k)/3-k(Bo./3)l, (15) 

p pp 

{uo.k ~Lj)f3 = (Bo.k)v-Bo(k.v). (16) 

These equations may now be manipulated to give a dispersion relation in the form ©(u, k) — 0, either directly or by noticing 
that they take the form of an eigenequation for eigenvalue cu, with eigenvector U""^ = (6'n, vj^, v^, /3^). 

Some factors of the dispersion relation can easily be found from equations (|l2|-p^. Taking the scalar product of k with 
(0), we find 

(c^-uo.k)/3.k = 0, (17) 

corresponding to an eigenvalue u = uo.k with eigenvector components /3 parallel to k and 9n = ~ = v = 0. This trivial 
eigenvector results from the requirement that V.B = be maintained, and has an amplitude identically zero. 
A further factor can be found, for which 6n ~ = 0. Using equations (p"2|-p^, we find that 

Bo.kf 



{iu - Uo.k) 



pp 



. ,2 (Bo.k)^ 

Uo.k) - ^ — 

pp 



PnLo{u; - Uo.k) } = 0, (18) 



with the components of the eigenvector having v parallel to Bq x k, 
iXp 
ijj + iXp 



(19) 



^ = -^^^v. (20) 
LU — Uo.k 

In these modes, two Alfven waves in the ionized gas are coupled to a shear mode in the neutral gas. We will refer to these as 
group A modes. 

The remaining 7 roots of the dispersion relation, which we will refer to as group B modes, are coupled by the remaining 2 
components of the drag force. To simplify the form of the dispersion relation somewhat, we assume (without loss of generality) 
that the wave vector k is parallel to x, and introduce scaled variables and parameters as follows 

"-p"'-T'"" = t'' (21) 

and a vector, va, with magnitude equal to the ionic Alfven speed and direction parallel to Bq, 
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Table 1. Properties of the most unstable solutions in Figure Q The panels of this figure corresponding to the various values of Ux are 
noted. Values are bracketed where the most unstable modes are not continuous with the resonant modes. 





Panel 


max[lm{uj)] 


kx at max[Im(a; 


0.22 


(b) 


0.020 


33.1 


0.43 


(c) 


(0.043) 


(2.70) 


0.78 


(d) 


0.516 


173.8 


1.22 


(e) 


(0.390) 


(6.08) 


1.43 


(f) 


(0.253) 


(4.06) 


2.43 


(g) 


0.163 


92.9 



0] 



VA = (22) 



Note that each of these variables has dimensions of velocity. This proves useful when interpreting the results of our analysis. 
In terms of these variables, the dispersion relation for the group A modes, equation (^8|), takes the form 

[(f2 — Ux)^ — VAx] ^ + i^i [(f^ — Ux)'^ — UAi] + ii^n^{i^ — Ux) = 0. (23) 

The dispersion relation for the group B modes, obtained from equations (p"2|-p^), has the form 
(fi - Ux) [(n - Ux){n + ifli) +imn] {{^'^ - cl + iQQi) [{fl-ux)^ - (fi + c?)] +ifin(^^ - Ux){if - cl)} 

-niQnclvAxl^VAx -VA.u]+cfvi^{n + ini){n'^ - cl+inQi) = 0. (24) 

It is useful to write this latter dispersion relation in the form 

Po{n) + i [niPi{n) + nnP2{n)] - nfPnin) - ninnPi2{n) - nlP22{n) = o, (25) 

where 

Po{n) = ~ c^) [{n ~ ux)^ - Vf] [{n ^ ux)^ - v!] , (26) 

Pi(n) = {2Q^ - cl) [(Q - ux)^ - vf] [{n - Ux)"" ~ v!] , (27) 
P2(n) = Q{n^-r^){Q-Ux)[2{Q-Ux)''-{vi + cf)], (28) 

Pii(n) = n[{n-ux)^-vf][{n-ux)'-v!]=Q[{n-ux)*-{vi + cf){n-ux)^ + cfvi,], (29) 
Pi2(n) = 2n\n^ux)^ -ivl + cf)n\n^ux)~ci[{n-ux)'' ~nvL + VA.uvAx], (30) 

P22(n) = Q{n-Ux)\il^ -cl), (31) 

and 

vl = + c? ± Vi^¥+W^^^x} (32) 

are the fast and slow ionic magnetosound wave speeds. 



3 NUMERICAL RESULTS 



In this section, we present numerical results for the roots of the dispersion relations ( p^ ) and (^^, concentrating in particular 
on circumstances where these roots correspond to physical instabilities. The roots were calculated as the eigenvalues of the 
complex matrix corresponding to equations (|l^-p^, using the routine ZGEEVX from LAPACK ( Anderson et al. 1999| ), and 



were verified by comparison to direct solutions to the polynomial dispersion relations. We show typical results in Figures ^ 
and ^ 

Where the coupling is weak (i.e. at high wavenumbers) and when the roots of the frozen system are well spaced, we find 
that the neutral acoustic modes are damped as tj = oj" — iXp/2, the neutral shear modes as uj = io° — iXp and the ionized 



fast, slow and Alfven modes all as lu = a;" — iXpn/2, in agreement with the results in Section 4.2 



Close to resonances, where two or more modes have similar phase velocities, additional modes need to be taken into 



account. The group A modes never become unstable, as confirmed by the analysis of Section 4.3. Instabilities are, however, 
found for group B modes. In Figure ^, we show the values of cj/k — uo.k for the group B modes, for one particular set of 
values of A, p, a, pn, Cn and B, and direction k. The plots in this figure show the variation of the roots as the slip velocity 
is varied, while the value of the wavenumber k is increased through the set of plots. Numerous graphs have been included 
to illustrate the full range of topologies in phase space, and to allow the development of the instabilities to be followed from 
linear through nonlinear order in A. 
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Figure 1. Variation of eigenvalues of group B with Ux, for various k. In this example, B = (1, 0, l)y/Ji, A = 1, p = 1, pn = 10, c? 
c2 = 1, k = (1,0, 0)A: and u = (1,0, -l)u:r. 
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Figure 1 - continued 
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Figure 2. Variation of eigenvalues of group B with fc, for various u^. As in Figure m, B = (1, 0, A = 1, p = 1, pn = 10, c\ = 0.1, 

cJ5 = 1, k = (1, 0, 0)A; and u = (1, 0, —\)ux- The velocities Ux are chosen to select various mode resonances. The circled points show the 
short-wavelength limiting values. 
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For large k (or short wavelength), instabilities (where Im(u;) > 0) appear at resonances between the neutral sound waves 
and the ionized fast- and slow-modes, and also (with a rather smaller growth rate) at resonances between one neutral shear 
mode and the ionized slow-modes. As k decreases, the morphology of the phase space evolves, and additional modes begin to 
influence the instabilities of the resonant system. In particular, the roots corresponding to the slow modes of the undamped 
system merge at ~ 25, leading to a cusped structure at Re(aj) — for smaller k. The apex of this cusp corresponds to 
solutions at Ux — > ±oo. This structure is particularly notable since, as k decreases, the roots at increasingly large Ux become 
unstable (cf. Section 4.7). Nevertheless, comparing the panels of Figure |^, it is apparent that the unstable modes which are 
seen can all be related to resonant instabilities for large k. 

In a physical flow, waves with the full range wavenumbers k vary independently (at least in the linear limit), for fixed A 
and u. In Figure ^ we show the variation of the roots for varying magnitude of the wave-vector, k, at fixed Ux (i.e. each of 
the points shown on these graphs appear simultaneously in a single physical system, as each corresponds to an independently- 
varying wave mode). As k decreases, each mode moves from its value at fe — > oo (shown circled), at which its phase velocity 
is determined by the uncoupled equations. As these examples were all chosen with two wave modes close to resonance, the 
resonant modes show a characteristic behaviour for large values of k, diverging from their positions at ^ oo in a direction 
essentially parallel to the imaginary axis. In each case, as fc — > two of the solutions have Im(tj) —> \{pi + pn) ~ 11, while 
the remaining five solutions have Im(tj) — » 0. For intermediate k, at least one root is on occasion unstable for each of the 
cases shown in Figure ^b-g). While each of these cases corresponds to a different mode resonance, and it is clear that these 
resonances have an important effect on the development of the roots with decreasing k, only three of the resonances lead 
directly to instabilities. 

The properties of the most unstable modes in each of these cases are given in Table |l| the cases in which the resonant 
modes lead directly to instability have most unstable modes at the largest wavenumbers. That the other cases are unstable is 
a result of the broad range of velocities away from resonance at which the three unstable mode couplings lead to instability, 
rather than resulting from the closer resonances which are present in these examples. Analytic criteria corresponding to these 
results will be determined in Section 



4.4 



We have shown results for each of the values of k which are present for particular physical parameters, for a particular 
direction of the wavevector. In the Section 4.1.2| , we present a geometrical argument which allows us to generalize these results 
for the full range of wavevector directions present in a physical flow. 



4 STABILITY ANALYSIS 

In this section, we will study the stability of the solutions to the dispersion relations (|2^) and ( |25| ) analytically. First, we 
develop general stability criteria for dispersion relations with the general form of equation (^), including in particular a 
discussion of the case in which two or more of the undamped wave modes have similar frequencies. We then apply these 
results, and find that isolated modes are stable, and that the group A modes remain so for all A. For the group B modes, 
however, we find that several two-mode resonances lead to instability. We also study higher resonances, and the stability of 
modes in the long- wavelength limit. 

Our results conflrm the importance of mode resonances for flow stability, as observed in the numerical results of the 
previous section. 



4.1 General considerations 

4.1.1 Perturbation theory 

We first consider the limit in which the coupling parameter A is small, i.e., that the damping parameters Qi and f^n in the 
dispersion relations are small. It is easy to see that Po{Q,) may be written 

j'o(j^)=n(" (33) 

a 

where ilL"' = Va are eigenvalues of the uncoupled system (i.e., when A = 0, and so Qi — = 0). We now use perturbation 
analysis, writing the frequency of each mode of equation (^5|) as a power-series expansion in the parameter A, which is assumed 
small: 

n„ = +n^^^ + ..., (34) 

where typically fii"' oc A". 

If there are no degenerate or near-degenerate roots in the uncoupled system, then the first-order corrections may be 
written 
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^(1) ^ , niPi(ni°') + »nP2(nL°') ^3^^ 

From this equation, it is apparent that for the flow to be stable, the roots of Po(f2) and QiPi{Q,) + QnP2{^) must interleave 
(cf Appendix |^; [Whitham 197^ . Since the roots of Po{0,) interleave with those of Pi{Q) and P2{^) individually, it is clear 



that, as > 0, this is indeed the case. (Note that this also demonstrates that 0,iPi{Q) + QnP2{^) has six real roots, 

which will be important in our discussion of long- wavelength modes below). 

However, if there are two strictly degenerate modes of the uncoupled system, with SI — \ then equation can 
no longer be applied and we have to take into account terms quadratic in A in our perturbation analysis. In this case, 
equations imply 

Pr{n^^^)^P2{n^^^)^0. (36) 

We then obtain the following quadratic equation for the first order correction to the degenerate mode, Q'^^ , 

i^i'^f n ^ " + * [^^Pli^'i'^) + f^nP2(f^L°')] fiL'^ - nfPn(r!L»)) - QlP22i^l^°^) - r!if7nPi2(nL°') = 0, (37) 

where Pl{0,) — dPi/dQ. This equation has two solutions, which corresponds to the breaking of the degeneracy of the modes 
in the presence of a small perturbation. Note that even though the expansion was to second order in quantities proportional 
to A, the correction is of first order. 

If we have near-degenerate roots, then Va ~ Vp for some set of roots {a,/3, . . .}: we will refer to this circumstance as a 
resonance. As a result, we have to consider the difference(s) Va — Vp, a ^ as further small parameter(s). For a two- wave 
resonance, we have 



n (nf - n(°)) + i (j7iP{(nW) + n^Pii^ii^^)) 



iniPl{n)^>)+innP2{nr)-^tPll{^^^')-^nP22{nr)-^i^nPl2{^^^^ = 0, (38) 

where the root of the quadratic with the smaller real part should be chosen for continuity with the non-resonant case. 
This equation can be written in the form 

Q2(J^L'')+jQi(t^L'') = 0, (39) 

with 

Q2(j7i^') = {iii'^r n ("^'" - + ' n - "^°') - ^'PM^^'^) - ^iP22{^L'^) - ^APi2{^L'^) (40) 

7,77^{a,/3} 7,7#a! 

and 

Qlin^^^) = (r!iP;(nL°)) + nnP^(nL°))) Q^^^ + niPl(nf ) + J7nP2(nf ). (41) 

We can now apply the Hermite-Biehler theorem (see Appendix ^), which tells us that the system described by the 
equation 

Qu(ni'^)+iQn-i(nL'') = (42) 

is stable (that is the roots of the equation have negative imaginary parts) when the roots of the polynomials Q„ and Qn-i 
axe real and interleave (and their leading terms are positive) , and that if this is not the case the system is unstable. In the 
case of equation ( |39[ ) we have polynomials of first and second degree, and can easily obtain conditions for stability. For higher 
resonances, the expressions are clumsy, so it is better to find instability conditions for particular cases rather than apply the 
general expression. 



4- 1-2 GeometTy of resonances 

In the preceding section, we presented a set of general criteria for the stability of fiows with internal damping parameters, 
and suggested that mode resonances may play an important role. We now find general conditions for the presence of mode 
resonances in two-fiuid MHD flows in three dimensions, by comparing the geometry of phase diagrams for the magnetosound 
modes of the ionized gas component and the sound and shear modes of the neutral component. 

In Figure ^(a) and (b), we show the phase speeds for the fast, slow and Alfven modes, as a polar plot of ui/k as a function 
of the angle between k and the magnetic field B, while in Figure ^(c) and (d) we show the same for sound waves and a shear 
mode in a frame in relative motion as a function of the angle between k and uq. 
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Figure 3. Phase velocity profiles, v = u)/k, as a function of the angle between k and B or uq, plotted on polar axes, (a), (b) Solid - fast 
and slow, dashed - Alfven, [the dotted curve - sound - is included for comparison but has no dynamical role here], c; > v\ for (a) and 
Ci < fA for (b); (c), (d) Solid - sound, dotted - shear, with velocity offset, (c) is for subsonic flow and (d) for supersonic flow. 



It is clear that the passing of the phase curve for sound waves in Figure |3|(d) through the origin corresponds to the 
Cerenkov condition for the emission of sound waves by an individual particle moving at uq. In a similar manner, resonance 
conditions between modes in the different phases can identified by looking for intersections between curves on suitably scaled 
and oriented versions of these plots (remembering to keep the origin at the same place in the ionized and the neutral plot). For 
example, resonances will always occur between the slow and shear modes, unless uo is parallel to B and uq > va- Likewise, 
resonances between the slow and neutral sound modes will be present in general once the relative motion is faster than the 
neutral sound speed, and for some choices of parameters once the slip velocity is greater than the difference between the 
neutral sound speed and the Alfven velocity. 

For finite damping, the resonance condition is weakened: the lines shown in Figure ^ may be thought of as blurred, 
although as a side-effect additional modes must also be considered in the stability analysis. 

We have now determined general stability criteria, suggested the importance of mode resonances in determining stability, 
and demonstrated that they will be an almost universal feature of MHD flows with drift velocities between components. 
However, depending on the form of the frictional force, only some of these resonances will result in flow instabilities. In the 
following subsections, we will apply these results to the dispersion relations given in Section ^ to determine quantitative 
stability criteria. 
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4.2 Non-resonant modes 



Using relation (|35[), first let us consider the stability of the modes of equation when there are no resonances in the 
uncoupled system, and Qi and fin are small. We find that all the modes are stable. In detail: 

Neutral shear modes For neutral shear modes, with fia"' — 0, 

fiL^> = -iQi. (43) 

Neutral sound waves For neutral sound waves, with fi?' = ±c„, 

(44) 

Magnetosound and Alfven modes For the fast and slow magnetosound waves, with fi^"' = w^; ± vi^s, and also for the 
(group A) Alfven modes, with ni"' = i: vax, 

nW=-i%. (45) 

4.3 Group A modes 

The stability of the group A solutions can easily be analysed for all A (not necessarily small) , by applying the Hermite-Biehler 
theorem to the dispersion relation, equation (^). This equation couples two Alfven modes and one neutral shear mode, and 
can be written 

Qa{n) + iQ2{n) ^ 0, (46) 
where 

Q3m=[iii-uxf -vL]n (47) 
and 

Qiin) = Qi [{n - - vl^] + nnfi(fi - u,). (48) 

It is clear that both Q2 and Qs have their full complement of zeroes, and that the zeroes of Q2 interleave with those of Q3 
(so long as there are no strict degeneracies and fin > 0). Hence, by the Hermite-Biehler theorem, the group A modes are 
always stable. 

4.4 Two- wave resonances 

We now consider the stability of the flow close to resonances between various pairs of group B modes. In order to determine 
whether instability occurs, we will follow the method of section |4.1.l| . 

4-4-^ Ionic magnetosound-neutral sound resonance 

First, we consider the case of a resonance between a neutral sound mode and an ionic magnetosound mode, so 

acn ~ + 7«f,s, (49) 

where a and 7 can independently take values +1 or —1. Applying the Hermite-Biehler theorem to the corresponding equations 
(^^-(^l|), one finds that instability occurs when the following inequalities hold: 

VAx{vf^^ - V^i) [-CiVsjSgn{vAx) + 'y{aCnVAx - VA-U)] < (50) 

and 

m , n \2 ^ 4:Vt,sivl,-vli)/vAx 2 

(fii+fin) > 7 T -, -(Ux + -yVl^s ~ CtCn) ■ (51) 

— CiVsjSgn(VAx) + jiaCnVAx — VA.U) 

In the limit of strict degeneracy, the r.h.s. of the second condition goes to zero, so in this limit the condition is no constraint. 
Note also that equation ( |5o| ) is the condition that the r.h.s. of equation (|^) is non-negative. Equation (|5l|) is in effect a 
minimum condition on A/fc for instability to occur (implying instability for large damping or long wavelength), for modes at 
a given distance from resonance. Note that these properties of conditions (^^ and (|^) also apply to conditions (^^ and dssl ) 
below. 

We can now consider particular cases of these general relations. 
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(i) acn ~ iti + iHf. 

So long as iiai; 7^ 0, equation niay be written 



< cf . (52) 



VA.U 

'yvi aCn 

L VAx 

If va -L u then 

7Qi;f < — . (53) 

Cn 

If Ci = as well, then 

7Q < 0, (54) 
which is possible only if 70 = —1, i.e. when the sound and fast-mode waves are oppositely directed. 

(ii) aCn ~ Ma: + fVs- 

Equation (^^ gives us 

VA.U" 



7*^s 



UAs 

If Va -L u then 



>cf. (55) 



javs > — (56) 

Cn 

which is possible only if 7a = 1, i.e. where the phase velocities of the sound and slow-mode waves are in the same direction. 
Unlike the previous case, we cannot take ci = in this relation, because in that case the remaining slow-mode wave will also 
have a similar frequency. The resulting higher resonance will be considered below. 

4-4-3 Ionic magnetosound-neutral shear mode resonance 

We now consider the case when Ux + 'iVf^^ ~ 0. Applying the Hermite-Biehler theorem to equations (p9|-|4l|) gives the instability 
conditions 

[uf_s -f 7?;Aa:VA.u] (Wf,s - vl^() < (57) 

and 

ton , n ^2 ^ 4uf,e(uF,e -t's,f)VwA , , .2 

(2ni + ^ln) > o ■ (u^+'yvi^s) ■ (58) 

Wf s + 7UAa:VA.U 



As before, we consider two particular cases 
(i) + -yvf ~ 0. 

In this case equation (|57[) gives us 



Vf < - fVAxVA-U (59) 

If VA -L u then the interacting modes are always stable, 
(ii) Uj: + 'yvs ~ 
Now we have from yl57\i 



Vb > — JVAxVA-U (60) 

SO that if VA -L u then this mode interaction will always lead to instability for sufficiently large A, so long as no other modes 
start to interact. 

4.5 Degenerate slow magnetosound viraves 

If we have Ci ~ (which is often the case in astrophysical applications), then the two slow magnetosound waves are resonant. 



i.e. Ux — Vs ^ Ux + Vs because Vs ~ 0. This resonance is of a different nature to those discussed in Section 4.4, as it is between 
two modes which propagate in the same phase. 

When Ci is strictly zero, two slow magnetosound waves become degenerate with common eigenvalue fl^^ = Ux. In order 



to 



investigate stability in this case, we have to solve equation (p7[), and find 



„(1) .fin , i /„, , -ifliQnClvAxivAxUx - VA-U) 

=''-^ 2 V + uxvliul-cl) ■ ^''^ 
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Instability can occur for the mode with positive choice of sign if 
vaxUx{uI - cl)[vAxUx - va-u] > 0, (62) 
For small but finite Ci, for instability to occur then in addition to condition (^2|), we also require that 

CnVAx[VAxUx - VA.uJ 

For the important case va -L u, condition ( |6^ leads to the simple relation 

ul>cl, (64) 

so long as vax 7^ 0. Note that the perturbation analysis will not apply when is close to c^, since Q^^ given by equation ( |6l| ) 
will be no longer be small. While no resonance with modes in the neutral phase has been assumed, this condition clearly 
relates to a resonance effect with the neutral sound waves in some fashion (indeed, when followed in numerical solutions, the 
mode is continuous with instabilities which originate at the slow magnetosound/neutral sound resonance for finite Ci). As the 
ionized sound speed becomes small, the system of coupled slow magnetosound modes becomes sensitive to the presence of 
other modes, even those rather far from apparent resonance conditions (cf. the sensitivity of the slow-mode waves to driving 



by other modes of nonlinear amplitude in the ionized gas, Falle & Hartquist 2002) 



4.6 Higher resonances 

We now consider the stability of some higher order mode couplings. 



4-6.1 Ionic slow magnetosound-neutral shear mode resonance 

Let us consider the case when Ux — ^ Q ^ Ux + Vs, which is possible if a, and Ux are small. 



First let us investigate the case c\ — Ux = 0, in which when we have a three times degenerate root Q}-^^ = for the 
uncoupled system. If we apply perturbation analysis to these roots, we obtain the following equation for the correction $1^': 

0. (65) 



'^A 

The polynomial Q2(fia ') has two identical roots flL^' = 0, which means that in order to satisfy the interleaving criterion, the 
polynomial Q3{Q.^a^) must also have one root — 0. This is possible only if va.u = - indeed for 



|vA.ul>2J^^ (66) 
V 27 VA 

Q3 has only a single real root. From this we conclude that the system is marginally stable when va.u = and unstable when 
VA.U 7^ 0. In the latter case, the equation has one term of second order, so the roots are 

^ ( "i"n.A.VA.u \ ^ ^g^^ 



\ / \ A / 

which scale as X^^'^ for A ^ 0, with one root being unstable. 

Now we consider the general resonant case where Ux and ci are small but finite. For the £7^°^ = mode we have 

ni'^ [{ni'^ - uxf - v!] + ^a [(nL^) - uxf - v!] + tn^i'^nl'^ - ux) - = o- (68) 

'^A 

The corrections derived for the roots which have fl^^ — Ux + jVa give equivalent results, relative to their alternative datum, 
as must be the case since the three roots of equation (^) correspond to the three resonant modes. 
The Q2 polynomial for equation ( |68| ) has two real roots 

(1) _ (20; + nn)ux ± ^/nlul: + 4»i(!:ii + nn)vS 

" " 2(0; + fin) • 

For the Qa polynomial, we first treat the case va-u = 0. Here, there are three real roots, one of which is zero and other 
two are 



ni''' =ux± x/vl + ninn^. (70) 



V 



Interleaving of these roots with those given in equation ( |69[ ) is violated (that is, the system is unstable) when 
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Vb < \Ux\ < , = — , (71) 

so we see that the system, which is marginally stable under fully degenerate conditions, can be unstable under resonant 
conditions. The range of unstable conditions widens in velocity as Qi and Sin increase. 
If va-u 7^ 0, we have an equation of third degree for Qa, 

- uxf ^ (v! + Q^Qn^] Qi'^ + = o. (72) 

Again the last term in this equation is of second order, so in the short- wavelength limit the roots are still given by equation (|67|). 
In the more general case where Q,iQn ^ VAi'VA-'n)^ /vXxy equation ( [t^ has complex roots (so the system is unstable) if 

{2ul - ISFi^i, + 27r2)'' >'i{ul+ ST if (73) 



(from the usual relations for cubic equations, e.g. Press et al. 1992), where we have introduced new variables 



ri = v! + Qinn^ (74) 



V 



r on VAxVA-U , , 

This implies that the solutions are unstable as va.u ±oo, other terms being equal: in particular, if is zero there are 
complex roots when 

|r2| >2(^)'^'. (76) 
If |r2| < condition (Q is satisfied for Ux in the range 

- sgn(r2)rl/^ - < ux < -sgn(r2)ri/^ + (^^^ . (77) 



If inequality ( |73[ ) is violated, then equation ( |72[ ) has three real roots and we have to solve the problem of interleaving. 
Criteria for interleaving may be obtained by inserting the roots of the Q2 polynomial, equation (^^, into the l.h.s. of 
equation ([t^). These are rather complex criteria to apply, except numerically. In general, however, for small Qi, fin the flow is 
unstable for conditions along whichever of the lines Ux = ±Vs satisfies equation (^o|); as Sli, fin increase this region of instability 
broadens and also expands along the other of Ux = ±Vs, in agreement with the limiting cases given by equations (^ij) and ([ttI). 

4-6.2 Neutral sound-neutral shear mode resonance 

Let us consider the case aCn ~ 0. One can show that equation for the corrections has the form 

nL'' W^^f - cl] + in, [2{ni''>f - d] - = o. (78) 

It is obvious that roots of Q3 and Q2 always interleave, and so the system is always stable under these resonant conditions 
(as is clear physically). 



4-6.3 Five-wave resonance 

There is also the possibility of a five- wave resonance if ~ aCn ^ Ux -\- 'yvs- But this case gives us a fifth order polynomial, 
which would be extremely complicated to analyze. So here we restrict ourselves to investigating only the case of full degeneracy, 
e.g. Cn = 0, Ci = and Ux — 0. Then we have three roots S7 = 0, and for the remaining two we have the equation for the 
corrections 

{ni^'>f + in'^^'>{2ni -\- ao - nii^i -f fin) = o. (79) 

It is clear from this that the system is stable. However, as we have seen above, the marginal stability of the three roots = 
under degenerate conditions does not guarantee stability under resonant conditions. 



4.7 Long-wavelength stability 

We now consider the stability of modes in long wave limit, i.e. for small wave-numbers. Note first that if we take the limit 
fc ^ in equation (^ ) then five modes have behaviour ui 0{k) and two modes have behaviour lo — > —iX{pi -\- pn) + 0{k) 
(cf. also the behaviour of the roots for long wavelengths shown in Figures ^(1) and 
Now we have to find the corrections to these limits of order k: 
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uj = Vck, a = 1...5 (80) 
and 

u = ~iX{pi + pn) + Vak, Q = 6, 7. (81) 



According to the definition of tlie limit, small k means that the second term in equation (|81|) is much smaller then the 
first one. As the first term has a negative imaginary part, the smaller second term cannot change the overall sign. Hence the 
roots given by equation (^l|) are stable in long wavelength limit. 

It is more difficult to determine the stability of the roots given by equation (^o|). To first order in k, the stability of these 
five modes in the limit of long wavelengths is determined by the roots of 

n-Pii{n) + nin,,Pi2{n) + nlP22{^) = o, (82) 

together with the condition that the leading term in this equation has a positive coefficient (which is the case). This equation 
is independent of k, and its roots are the velocities Va in equation (|8o|). As equation ( ^^ is a polynomial equation with real 
coefficients, complex roots can arise as conjugate pairs and any such complex roots in equation (^2|) will result in instability 
in the limit of small k. 

Real roots of equation (^^ are marginally stable to first order. To determine their asymptotic stability, we would need to 
study the interleaving of these roots with those of QiPi{Q) + QnP2{^)'- if the roots of these two polynomials do not interleave, 
the asymptotic solution will be unstable at order k^. Note that we know that fliPi{Q) + Q,nP2{^) has its full complement of 
real roots, from the discussion after equation ( ^5| ) above. 

We will first study the stability of the long-wavelength solutions in various limiting cases. For u = we have five real 
(and hence stable) roots, which correspond to the fast- and slow-mode waves of the fully coupled system together with an 
additional mode with f2 = 0, as would be expected on physical grounds. 

If Ux is large compared to any other characteristic velocity, the solutions are f2 = [or more accurately Q,nC^/ {^iUx)], 
and QiUx/iJli + 0,n) and Ux (the latter two are double roots). This can be easily seen if we write equation ( [s^ as follows 

- Uxf [(fii + fin)!^ - ^iUxf -iSl- Ux) [(!^i -I- 0,^)0. - 0.iUx\ [^i{vi + C?,)Q. + Q.r,cl{Q. - Ux)] + 

^^,o,v\^^l + Q,;Q,^c\vp^x[vAx^-^A-'a\ = 0. (83) 

The presence of double roots means these solutions are only marginally stable in the limit, so we need to carry the analysis 
to higher order to determine stability in the asymptotic regime. To lift these degeneracies, we should expand to lower order 
in fl and Ux, or to higher order in k (or include nonlinear terms neglected in our initial linearization of the problem). 

If k _L Va, i.e. vax = with all other terms finite, equation ( p3[ ) has roots Q — Ux and — QiUx/(Oi + 0,n)- The remaining 
equation of three roots satisfy 

mn [{n - Uxf - (vl + c?)] + Q„(n - ux){n^ ~ cl) = o, (84) 

where ai = +Qn), Qn ~ 0,n/{^li +Q,n)- Considering the ordering of the roots of the factor proportional to ai and that 

proportional to Qn = 1 — cvi, it is clear that there must be at least two real roots of equation (^^ for all < Qi < 1, and hence 
there will be three real roots. Therefore for vax = 0, all the solutions are stable to order 0{k). 

If we consider corrections to the 0{ux) asymptotic solutions, then for the two roots Q — Ux + S, we find that the highest 
order corrections have S = 0{l/ux), and are given by 



QnCivAx'VA-U/Ux — v\^(QnCl + Qicf) 



(85) 



which can clearly be complex, e.g. if va.u — and ci = the criterion for instability becomes ^vj^^c^/vX > Sli/fin. Similar 
expressions for the case ~ QiUx/{Qi + Qn) are more complex, but the condition for instability is 

4f2f [wi^r2i(r2icf + Qncl) - Qni^i + 0.n)clvAxVA-U/Ux] > [0,f {v% + cf) + ^.^clf . (86) 

If, instead, we include higher order terms in k, where in the limit of large Ux 

Po = Q^{n~ux)\ (87) 

niPi + QnP2 = 2Q'^{n-ux)^[{ni + nn)^-uxni], (88) 
nf Pii + ninnPi2 + ftlP22 = ^^(^7-u^)^[(fii + nn)n-^t,^^i]^ (89) 

then the common roots S7 = and Ux (twice) are marginally stable. Interleaving is satisfied for the remaining four roots, so 
they are each stable for all k at large Ux (although only the limiting form as fe — ^ is required for the present discussion). 

First-order instabilities can also arise in the intermediate velocity regime. Looking again at equation (^), we have shown 
that quintic consisting of the first two terms has five real roots. However, the values vax and va.u can be varied independently 
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Figure 4. Values of H for limit of long wavelengths, for da = 1000 km s~^. (a) and (b) show the real and imaginary parts of two roots 
which are complex for part of the velocity range. The remaining three roots are stable and not shown here. The roots do not tend to 
their asymptotic behaviour in this plot, as the maximal slip velocity Mmax = 100 is small relative to the Alfven velocity (cf. the case of 
vx = 20kms~^ in Figure Bl). 



of any of the other parameters of the overall equation, so the additional linear function consisting of the third and fourth 
terms is entirely arbitrary. Even when va.u = 0, the number of real roots of equation ( ^ ) can be as few as one. 

As the long wave limit is quite important for the astrophysics of molecular clouds and stellar winds, in the next section 
we will study a typical astrophysical example, SiO maser spots in late-type stellar winds. 



5 SIO MASER SPOTS IN LATE-TYPE STELLAR WINDS 



As mentioned in the introduction, the current result s are applicable to mass loss from highly evolved stars, which are believed 
to be driven by radiation pressure on dust particles ( Bowen 1988 ; MacGregor fc Stencel 1992 ; Mastrodemos, Morris fc Castor 
1996 ; 3imis, Icke fc Dominik 2001 ). The winds are subject to strong perturbations as a result of stellar pulsation, and models 
have found both radiation-pressure instabilities and others resulting from the condensation of the dust grains, as well as 
radiative instabilities. The structures which result from these instabilities may be imp ortan t in determining the properties of 
planetary nebula halos which will form as the star ages further. Mastrodemos et al. ( 1996 ) find that slip between dust and 
gas components leads to fo rmation of dense shells at intermediate radii, but that these dissipate as the wind moves away from 
the star. Simis et al. ( 2001 ), however, find that dense shells survive to large radii, and that an accurate treatment of dust-gas 
slip is essential in treatments of late-type winds. 

If, in the wind of a late-type star, ion-neutral friction is the main source of coupling and the streaming speed is less than 
about 30kms~^, A can be calculated from the results obtained by Osterbrock (L961) for the ion-neutral momentum transfer 
cross section, so 



A ^ 1 X 10- 



( —!—\ ( 

UO-s^ I Pi 



(90) 



where x is the fractional ionization, i.e. the ratio of the number density of ions to the number density of hydrogen nuclei. 
The value of k in the dimensionless units used in Figures ^ and ^ is related to the wavenumber in physical units fcphys by 
k = fcphys/fco, where 



fco 



%/2Api 



3 X 10-^^ cm-^ (loii'lVa) (io3kms-) ( iF^) 



(91) 



VAi Vl0"cm-3y VlO^kms- 

To compare this with the scale-lengths of features in the winds, we must have some idea of the appropriate values of 
VAi, Un, X and of the scale-lengths of interest. SiO masers are found in the outflows of some evolved stars at positions near 
those at which dust is expected to form (e.g. Kemball & Diamond 1997 ). The magnet ic field strength in the maser spots 



can be inferred with difficulty from the observed pola rization (e.g. Watson fc Wyld 199^ ), although it is not certain that the 
polarization is in fact caused by the Ze eman effect (VViebe fc Watson 199S). T he resulting estimates of the magnetic field 
strengths are in the range of 2 to 10 G (lElitzur 1996| ; [Kemball fc Diamond 1997]). Pumping models for these masers suggest 



that the neutral density in the spots is n-n ~ 10"^ - 10" cm ([Doel et al. 19951) . If one percent of the mass is contained in 



grains, vai could be as high as 10"^ kms ^, but may be closer to 100 km s ^ . The size of individual maser spots is ~ 10^^ cm 
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Figure 5. Values of £7 in tlie long wavelength limit, for vp^ = 20kms~^. (a) shows the real parts of the roots for small Ux- The regions in 
(a) where the real parts are distinct correspond to stable roots with zero imaginary parts: (b) shows the imaginary parts of the complex 
conjugate pairs, (c) and (d) are the same as (a) and (b), except that they are plotted for a rather wider range of Ux, and the asymptotic 
behaviours of the roots are shown as dashed curves. At large Ux, the roots with the largest real parts follow the expected asymptotic 
behaviour, but have finite imaginary parts. The next two roots should converge to a common asymptotic limit, Q ~ O.OlUx, but while 
the variation of the smaller of this pair is almost indistinguishable from this form (indeed, it overlies the dashed curve showing the limit 
on the graph), the larger lies well away from its asymptotic value, and convergence is only clear at far larger Ux, as shown in (e). 
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(e.g. ^emball fc Diamond 1997| ), for which fc « 3 x 10"" cm if the maser spot size corresponded to a half wavelength of a 



perturbation. The fractional ionization is uncertain, but it is unlikely that grains will form where carbon is ionized, and some 
of the elements with low ionization potentials will be depleted substantially from the gas phase due to the grain formation 
process. Grain-neutral friction for a grain-neutral relative speed comparable to the thermal speed of the neutral material is 
of similar magnitude to the ion-neutral friction for x between lO"'^ and 10~^, if the grains have a fractional abundance and 



size distribution similar to those in the interstellar medium (e.g. Baker 1979). 

Therefore the SiO maser regions are probably large enough that their overall properties correspond to our long wavelength 
limit. To study their stability in detail, we consider parameters characteristic of the regions of late-type stellar winds with 
SiO maser spots, as follows: pi/pn ~ 0.01, Cn — 3kms~^, ci = 0, ua = 1000 kms"'^ and u = 0-100 kms~^. We also assume 
that u is perpendicular to magnetic field, i.e. va.u = 0, and in particular study the case where the wave vector, k, bisects the 
angle between u and va. Numerical solutions of the long- wavelength dispersion relation are shown in Fig. ^. For comparison 
we also show in Fig. ^the calculations for smaller Alfven speed = 20kms~^. Note that for some slip speeds there are two 
pairs of complex roots, but that for larger Alfven velocities the higher velocity roots are stable at all slip speeds. We will refer 
to the instability of the higher velocity roots, shown in Figs ^,b, as a type I instability, and that in Figs ^ and ^,d as type 
II. 

It will be noticed that as ua changes from 20 to 1000 kms~^, the region of type II instability moves to higher slip velocities. 
For ua ~ 20kms~^, the maximum of Im(n) is around Ux ~ 1.5kms~^ and for ha ~ 1000 km s~^ this maximum is around 
Ua; ~ 20 — 25kms~^. The value of the maximum growth rate, however, changes little, and is around 0.06-0.1 kms~^. 

For small Alfven velocities (like 20 kms"^), a type I instability is also present. The regions of two instabilities overlap 
to a degree, with the range of type II instability being for — 1.2-1.6 km s"^ and that for the type-1 instability being for 
Ux <; 1.5 kms~^. 

The maser spots are typically at distances of 5-10 a.u. from the centres of the stars. For an outflow speed of the order of 
lOkms"^, we conclude that an instability will have no significant effect on the maser spots unless it grows on timescale of a 
few years or less. From Fig. ^b one can see that for relatively small Alfven velocities this is indeed the case, and such maser 
spots would be destroyed by type I instability. But as maser spots exist at these radii, we can conclude that magnetic field 
(and therefore the Alfven velocity) cannot to be too small. Indeed, for Alfven velocities around 1000 kms~^, type-I instability 
is suppressed and we only have to deal with type-ll instability (see Fig. This instability at its maximum grows on timescale 
of order 5-10 years, but the range of streaming velocities at which the instability has even this large a growth rate is quite 
narrow, between 28 and 35kms~^. The streaming speed of grains through a maser spot is uncertain and depends on the 
fractional ionization and the grain size distribution; it may be that in reality the streaming speed is too large or too small 
for the instability to be driven at a rate high enough to destroy the maser spot. It will be important, as models of late- type 



stellar wind including both the chemistry of grain formation and dynamics are developed (e.g. MacGregor & Stencel 1992 



Cherchneff 1995; Gail & Sedlmayr 1998), to determine the relationship between the calculated grain-neutral streaming speed 
and the input magnetic field. These results are necessary for a full application to stellar outflows of the present analysis. 

Finally, it is of interest to determine whether the short wavelength analysis developed in previous sections can be applied 
to the flner-scale properties of the maser spots. This will be the case if the widths of the resonant interactions are narrower 
than the spacing of non-resonant modes, i.e. that 

^ I4>in, (92) 

where Knin is the smallest characteristic velocity of the problem which should be ~ lkms~^, and we take pn in the left 
hand side as it is larger than pi. This inequality means that, in this example, short wavelengths are 27r/fc ^ 10^" cm. 

This maximum wavelength must be compared to minimum wavelengths required for the fluid approximations to be valid, 
and for the growth of the two-fluid instability to outweigh wave damping as a result of viscosity. The maximum mean free 
paths for both electrons and ions are 



Ac i ~ 10^ cm X 



r 



) (lo^^)" 

(Melrose 1986), while for neutral atomic species, the free path is 

An ~ — (94) 

- ^x^°'^-x(lO^^)"' 

assuming a geometric cross-section. For longer wavelengths, the finite collision frequencies will lead to viscous damping of waves 
in a number of periods given by the ratio between the wavelength and the free path, and any instability must grow rapidly 
enough to outweigh this damping. For a given slip velocity, the growth rates for the instability we discuss are proportional 
to the density while the viscous damping rates are proportional to the ratio of the square of the frequency divided by the 
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density, so damping will indeed only be important for high-frequency waves in diffuse media. In fact, at these smaller scales 
further instability modes may also become important (driven by, e.g., the Hall effect, Balbus & Terquem 2001). 

For the particular example under discussion, we require 10^ cm ^ Aw ^ 10^"^ cm which obviously can be satisfied, so 
that our short wavelength analysis is applicable for a substantial range of wavelengths. In general, an appreciable range 
of potentially unstable wavelengths remains for all temperatures T <C 10^ K. While the shorter wavelengths have less direct 
observational significance, their growth to nonlinear amplitudes will likely have substantial effects on the overall flow structure. 



6 DISCUSSION AND CONCLUSIONS 

We have seen that the slip between different fluids in a multifluid medium can drive instabilities. At short wavelengths, several 
modes of instability result from resonances between waves propagating in the different fluids; these instabilities remain for 
longer wavelengths but they become increasingly inter-coupled and less easy to characterise. In particular, neutral shear/ionized 
slow-mode resonance instabilities will grow in almost all cases where there are finite drift velocities. The maximum growth 
rates we find are of order the characteristic ion/neutral collision rates in gas without internal slip, which are typically 

:/in = 1.3 X 10"^°nns~^ (96) 



(Osterbrock 1961; Brandenburg & Zweibel 1995), rather than inversely proportional to density as is the case for instabilities 
derived from leakage of magnetic fiux as a result of ambipolar diffusion. This corresponds to a characteristic lengthscale of 
2.5 X 10^^(11/1 kms^^)7in^ pc, where we note that the most unstable wavelengths may be one hundr edth of this (Table p' ). 



The instabilities appear to be similar to the two-stream instability of plasma flows (e.g. Melrose 1986; Bingham et al. 2000). 

In Section ^, we apply our analysis in its long-wavelength limit to the overall properties of SiO maser spots in late-type 
stellar winds. We flnd that if the SiO spots are not to be subject to violent instabilities, the slip velocity between the phases 
in these regions must lie outside certain limits, or the magnetic field in these spots must be very strong. 

Our short wavelengths results also have important implications for many other astrophysical systems. In molecular 
clouds, we expect that the ionized sound speed is small and that the magnetic field is perpendicular to the inter-component 
slip velocity, but that this slip velocity is not small. This means that we always can flnd directions such that inequality ( [7l[ ) 
is satisfled. It seems likely, therefore, that molecular clouds are generically unstable to the growth of slow-mode waves. In 
the present analysis we can only conjecture the non-linear endpoint of these instabilities, but the two obvious possibilities - 
fractionation of the phases as a 'slugged' flow, with consequent rapid loss of magnetic field support, or the limiting of the wave 
spectrum at finite amplitude - each have clear practical and observational consequences for the ecology of the interstellar 
medium. 

Our present study is substantially simplified. This has allowed us to derive some rather general results. In future work, 
we will model the structure of the systems of interest more completely, including the spatial structure of the background 
flow, more interacting phases (e.g., treating electrons, ions, neutrals and a spectrum of sizes of dust particle as independent 
species), variation of the frictional constants with state, and the detailed nonlinear evolution of the instabilities. 
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APPENDIX A: THE HERMITE-BIEHLER THEOREM 

We quote here witho ut pr oof the Hermite-Biehler theorem, a very general result which underlies much of our analysis. From 
theorem 4' of Levin ( |l964| ), Chapter VIL 

Theorem 1. (Hermite-Biehler) In order that the polynomial 

w{z) = u{z) + iv[z), (Al) 

where u{z) and v{z) are real polynomials, does not have any roots in the closed upper half-plane \xa{z) > 0, it is necessary 
and sufficient that the following conditions are satisfied 

(i) the polynomials u{z) and v{z) have only simple real roots, and these roots separate one another, i.e., between two 
successive roots of one of these polynomials there lies exactly one root of the other; 

(ii) at some point xo of the real axis 

v'{xo)u{xo) ~ v{xo)u'{xo) < 0, (A2) 
where u'{z) = du{z)/dz and v'{z) = dv{z)/dz. 
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